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The network formed by ridges in a straightened sheet of crumpled paper is studied using a laser 
profilometer. Square sheets of paper were crumpled into balls, unfolded and their height profile 
measured. From these profiles the imposed ridges were extracted as networks. Nodes were defined 
as intersections between ridges, and links as the various ridges connecting the nodes. Many network 
and spatial properties have been investigated. The tail of the ridge length distribution was found 
to follow a power-law whereas the shorter ridges followed a log-normal distribution. The degree 
distribution was found to have an exponentially decaying tail, and the degree correlation was found 
to be disassortative. The facets created by the ridges and the Voronoi diagram formed by the nodes 
have also been investigated. 

PACS numbers: 89.75.Hc,83.60.-a,89.75.Fb,42.62.-b 



INTRODUCTION 

The crumpling of paper is an everyday occurrence, yet 
it is a surprisingly rich and complex process. Paper is 
an elastic, flexible and heterogeneous material and many 
authors have tried to describe its crumpling properties 
analytically El B, §|, numerically 0, 0| and experimen- 
tally [rj M.MM]iOTl11|. The crumpling process of paper 
is also interesting because it is a special case of the thin 
plate deformation problem that is central in describing 
processes that occur for example in car crashes and tank 
failures 12j. Earlier studies have tried to describe the 
ridge network of crumpled paper 0, @, E3] and some re- 
sults have been found, however much is still unclear. This 
work aims at describing the ridge network formed dur- 
ing a common hand crumpling process of ordinary print- 
ing pap er. The application of modern network theory 
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151 ] have been specially emphasized. 
This paper is organized as follows. In section 2 the 
experimental procedure is described, and in section 3 the 
ridge detection method is presented. The results are dis- 
cussed in section 4, in particular the ridge length and the 
degree distribution are discussed. Also the degree-degree 
correlation, the clustering and the surface roughness is 
investigated in addition to the facet distribution and the 
angular ridge distribution. Finally the main conclusions 
are summarized in section 5. 



EXPERIMENTAL PROCEDURE 

Ordinary printing paper was used for all the experi- 
ments, and some of the properties of the paper are given 
in table [J All the samples were cut into square sheets of 
21 x 21 cm, and crumpled by hand into small balls. The 
diameter of the various balls produced are given in ta- 
ble HI The hand crumpling procedure have been applied 



before [fj, la, |lfjj, and is practical because it is easy to 
conduct and produces a compact result. Unfortunately 
the process in not repeatable and poorly controlled. Sev- 
eral test crumplings were conducted before the measured 
samples were crumpled in order to reduce the variance 
between the samples. Earlier studies [8] on acoustic emis- 
sions from crumpling of various materials have indicated 
that the emission spectra show a surprisingly low sensi- 
tivity to the crumpling method. This may indicate that 
the outcome of the crumpling is not highly sensitive to 
the details of the process. Balankin et al. [1(3] discuss 
the scaling behavior of the crumpling process for differ- 
ent paper thicknesses. They conclude that the impact 
of the variation of the applied confinement force F on 
the ball radius R is small since there is only a weak 
dependence R cx F~°- 25 , For these reasons no special 
precautions, such as dents or initial folding, were taken 
to increase repeatability. After crumpling the samples, 
they were carefully unfolded taking care not to tear the 
paper, introduce new ridges or remove some of the origi- 
nal ridges. When the paper ball was unfolded the paper 
was stretched to a size of 20 x 20 cm, and fastened to an 
aluminum plate. This ensured that the vertical height of 
the samples were no more than 12 mm (the maximum 
range for the instrument used). 

The full (2+l)-dimensional height mapping was mea- 
sured profile by profile using a laser profilometer over an 
area of 18 x 18 cm in the center of the samples. The 
height of each point was measured using a laser giving 
a voltage output linearly proportional to the distance 
between the probe and the paper surface. The voltage 
output was converted to a floating point length-measure 
using a 16 bit AD converter. The laser diameter used was 
30 /im, however accuracy considerably smaller than this 
could be achieved. Each profile was acquired by sliding 
the sample under the probe while measuring. Multiple 
profiles were acquired by stepping the probe normal to 
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Sample X step Y step Thickness Weight Ball diam. 









[ fim] 


[g/m 2 ] 


[mm] 


1 


900 


900 


51 ±5 


49.0 ± 1.0 


26 ± 2 


2 


1800 


1800 


51 ± 5 


50.0 ± 1.0 


27 ± 2 


3 


1000 


1000 


95 ± 5 


80.0 ± 0.5 


35 ± 2 


4 


900 


900 


95 ± 5 


80.0 ± 0.5 


32 ± 2 


5 


3600 


3600 


100 ± 2 


83.0 ± 0.5 


33 ± 2 


6 


900 


900 


220 ± 5 


175.0 ± 1.0 


43 ± 2 



TABLE I: List of samples investigated. X steps and Y steps 
is the number of measured points in the X and Y direction 
respectively. Thickness is the thickness of the paper and ball 
diameter is the diameter of the ball produced during the crum- 
pling process. All samples was originally 21 x 21 cm and 
thereafter unfolded and stretched to 20 x 20 cm producing a 
maximum height of 12 mm. An area of 18 x 18 cm in the 
center of the samples were measured. 



Sample Nodes Links C Cd Cr Max Deg. 



1 


503 


890 


0.182 


0.4371 


0.0045 


8 


2 


1211 


2238 


0.190 


0.4315 


0.0020 


9 


3 


190 


293 


0.138 


0.4458 


0.0095 


6 


4 


350 


580 


0.162 


0.4394 


0.0064 


8 


5 


929 


1829 


0.231 


0.4326 


0.0029 


10 


6 


286 


501 


0.199 


0.4384 


0.0083 


8 



TABLE II: List of extracted networks with their number of 
nodes, number of links, clustering coefficient, C, the cluster- 
ing coefficient for the corresponding planar Delaunay network, 
Cd, the clustering coefficient for a non- planar randomized 
network with the same degree distribution, Cr, and the max- 
imum node degree for the network, Max Deg. The clustering 
coefficients for the random networks was calculated using an 
average over 1000 samples after each sample had 10000 ran- 
dom rewiring. 



the sliding direction. A typical one-dimensional height 
profile and a complete (2+l)-dimensional map are shown 
in figure [T] The number of points per profile was kept 
equal to the number of profiles resulting in a square grid 
of measurements or "pixels" . The number of points used 
for the various samples are given in table |TJ The in-plane 
accuracy of each point was no larger than 10 /im for any 
sample, and in the out of plane direction it was 0.5 /im 
for all samples. 

RIDGE DETECTION 

A ridge stands out as a line of high curvature in a 
otherwise smooth landscape. The curvature of any point 
in a height profile can be calculated as the field 

V 2 £(aO, where x — (x,y) is the planar coordinates. In 
the present case it proved necessary to smooth the height 
profile with a short range gaussian filter before cal- 
culating \/ 2 £,{x) in order to filter away small scale fea- 
tures. A range of different filters were tested, and the 
result did not seem sensitive to the details of the filter. 
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FIG. 1: (Color On-line) Top: Grey scale height plot of sam- 
ple 6 as function of x and y position. Color indicate elevation 
(brighter is higher), ridges are clearly visible. Middle: A sin- 
gle one dimensional profile from sample 6 (marked as a black 
line in the bottom plot). The points of the profile that give 
rise to ridges are marked by squares. Bottom: Network ex- 
tracted from sample 6 superimposed on the gray scale plot 
from the top figure. 
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The main effect of the filtering was the removal of single 
isolated high-curvature pixels, or small groups of such, 
and a narrowing of the ridge lines. After filtering, the 
curvature field was calculated and thresholded so that 
all points over a given value were considered to have a 
unit value and all other points to have a zero value. Any 
isolated points above the threshold were filtered away. 
From the remaining points lines were detected as ridges. 
It is throughout this paper assumed that all ridges are 
straight lines. It proved difficult to automate the ridge 
extraction process from the thresholded field, finally this 
step had to be done manually. Some statistics of the pro- 
duced networks are listed in table HH Figure Q] shows an 
example of a full ridge network. In the middle plot of 
figure Q] a single one-dimensional profile is given, and all 
points along this profile giving rise to ridges are marked. 
It can be noted from this figure that not all sections of 
the profile that have high curvature give rise to a ridge, 
while some smooth sections does give rise to a ridge. This 
may stem from the directionality of the ridges relative to 
the profile shown. Ridges crossing the profile at a small 
angle may seem smooth, but small local dents crossing 
close to orthogonally may seem large. 

Nodes are defined as intersections between ridges, and 
a ridge therefore only extends from one node to another. 
All the links are regarded as undirected since a paper 
ridge does not have any preferred direction. The net- 
works formed are fully connected and have therefore only 
one component. 

NETWORK PROPERTIES 

The different paper thicknesses used in the experiments 
showed a clear trend that thinner paper crumple more 
than thick paper, and therefore produce more nodes and 
links (see table U and |ll| . Apart from the scale of the 
network created, no significant differences in the various 
distributions referred to below was detectable. As a con- 
sequence most distributions are averaged over all sam- 
ples after each of them have been normalized appropri- 
ately. The lack of change in the behavior due to sample 
thickness may arise from the small amount of data avail- 
able, and no correspondence between paper thickness and 
other properties can be excluded. From a scaling point 
of view a qualitative change of behavior is not expected 
since a large and thick sheet of paper is equivalent to a 
thin and small sheet. Note that the needed confinement 
also varies with the paper thickness, and all our exper- 
iments are conducted at approximately the same con- 
finement. Sultan and Boudaoud [ll[ discuss two regimes 
for the crumpling process depending on the confinement 
of the sample. The transition confinement is partly de- 
pendent on the paper thickness. Our experiments are 
as mentioned conducted at approximately constant con- 
finement (although it is poorly controlled), and it might 



therefore be that due to the varying paper thickness our 
samples lie in different regimes. However the uniform be- 
havior of the samples indicates that they are all in the 
same regime. Also the number of self-contacts are very 
large for all the samples, and this indicates that they are 
all in the highly confined regime. 

It can be seen from table U that samples 3 and 4 both 
have the same paper thickness, although they have a sig- 
nificantly different number of links and nodes. This is 
most likely due to the difference in confinement. Sample 
3 had a larger ball radius than sample 4, and was there- 
fore less confined, and have also fewer nodes and links 
than sample 4. 

Ridge length 

The length of a ridge between node a and b is de- 
fined as the spatial length from node a to b, following the 
assumption that all ridges are straight lines. Previous 
works have reported log-normal, gamma and exponential 
functions 0, d, [l(| to give good fits for this distribution. 
However we find that, whereas the small scale part of 
the distribution is well fitted by a log-normal function, 
the tail of the distribution is not well fitted by any of 
the above mentioned functions. The large scale part of 
the distribution is better fitted by a power-law function 
p(l) oc (1 — l/loY where I is the ridge length and Iq is 
the maximum ridge length for a given sample. Both fits 
can be seen in figure [2l We have found the tail to be 
best fitted by an exponent (3 — 0.81. To compare the 
fits of the different functions they are plotted in figure [3] 
divided by the original distribution in order to emphasize 
any discrepancies. 

The underlying reason for the shift in behavior may 
stem from the fact that the distribution of short ridges is 
dominated by remnants of originally long ridges. These 
ridges have been intersected by " younger" ridges crossing 
them after their formation. As outlined by Blair and 
Kudrolli @ this random sectioning of ridges will give rise 
to a log-normal length distribution. The larger ridges 
on the other hand have not been so heavily sectioned by 
younger ridges. They are therefore not expected to follow 
the log-normal distribution of the shorter ridges. Instead 
we detect a power-law dependency of the distribution of 
the difference between the longest ridge Iq and the ridge 
length. It is reasonable to assume that larger samples 
will produce larger maximum ridges, and therefore lo is 
a sample size dependent quantity. Why this difference 
should exhibit a scale-free behavior is not clear. 

Degree distribution 

The degree of a node is defined as the number of ridges 
meeting at that node. The distribution has been found to 
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FIG. 2: (Color On-line) a) Plot of the average non-cumulative ridge length distribution, p(l), as a function of 1 — l/lo where 
I is the ridge length and io is the maximum ridge length for any given sample. The data is fitted by a log-normal and a 
power- law p(l) oc (1 - l/lo) with (3 = 0.81. b) Plot of the average non-cumulative ridge length distribution, p(l), as a function 
of normalized ridge length l/lo- 
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FIG. 3: (Color On-line) For comparison the results shown in 
figure [5] from fitting the ridge length distribution, p(l), with 
a gamma, log-normal and exponential function divided by 
the data itself are shown together with the same plot for the 
power-law fit as a function of normalized ridge length l/lo- 



FIG. 4: (Color On-line) Plot of the degree distribution, p(k), 
as a function of node degree k with a fitted log-normal tail. 
The insert shows the same data plotted on log-log scale. This 
shows that the distribution does not have a power-law dis- 
tributed tail. 



have a maximum probability at a median degree and pro- 
duce a gaussian like form that is plotted in figure 0J The 
tail of the distribution is well fitted by a log-normal func- 
tion of the same form as in equation 2J This is in strong 
contrast to many naturally occurring networks that show 
a power-law tail giving a larger portion of high degree 
nodes than can be seen in the acquired samples. 



Degree-degree correlation 

The correlation between the degree of connected nodes 
has been studied using the procedure developed by 
Maslov and Sneppen |15j. They have defined a corre- 



lation measure 



C{k u k 2 ) = 



P(k u k 2 ) 

P R (h,k 2 ) 



(1) 



where P{ki, k 2 ) is the probability that a node of degree 
k\ is linked to a node of degree k 2 . Pr(A;i, k 2 ) is the same 
average probability for a set of randomized networks. 
The randomized networks are assumed to have the same 
number of nodes and links, and the same degree distri- 
bution as the original network. A value C(ki,k 2 ) > 1 
indicates that there is an over-representation of links be- 
tween nodes with degree fci and k 2 , whereas C{k\, k 2 ) < 1 
indicates an under-representation. In order to look at 
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FIG. 5: (Color On-line) Plot of the correlation matrix 
C(ki,k2) = P{ki,k2)/Pii(k\,k2) for each sample. The plots 
indicate that the networks are disassortative. 



FIG. 6: (Color On-line) Plot of the Z-matrix Z(ki,k 2 ) = 
(P(ki, k2) — Pi?(fci, fc2))/<r«(A;i, fei) for each sample. As in the 
C(ki,k2) case of figure [5] the plot indicates a disassortative 
trend. 



the statistical signification of the correlation, Maslov and 
Sneppen introduced another correlation measure 



Z(ki, k 2 ) 



PjkiM -P fl (fci,fc 2 ) 

GR{kl,k%) 



(2) 



where an{k\, k 2 ) is the standard deviation of the sam- 
ples used to generate Pn(ki, k 2 ). For P(ki,k 2 ) only the 
sample data is available. If a given coupling P(ki,k 2 ) 
is over-represented (that is P(ki,k 2 ) > Pn(ki, k 2 )) then 
Z(k\, k 2 ) > and if it is under-represented Z(k\ 1 k 2 ) < 0. 
If the standard deviation is small the corresponding cor- 
relation coefficients are large, thus emphasizing statisti- 
cally significant results. In all results presented here 1000 
randomized versions of the various samples were used to 
produce PR(kx,k 2 ) and o~n{ki,k 2 ). Each randomization 
used 10000 rewirings of the original network. 

Figure [5] shows the C(ki,k 2 ) matrix for all the sam- 
ples. There is a tendency of small degree nodes not to 
link to other small degree nodes, but rather link to large 
degree nodes. Links between large degree nodes is also 
under-represented. This type of networks is known as 
disassortative networks. Figure [6] shows the Z(ki,k 2 ) 
matrices for the same samples, and the same trends as 
in figure [5] can be observed. There is a clear trend in 
nearly all examined networks [l6| that technical and bio- 
logical networks such as the Internet and various protein 
interaction networks are disassortative, and that social 
networks such as acquaintance networks are assortative. 
The underlying reason for this is still not fully under- 
stood. 



Clustering 

The cluster coefficient C for all the samples are given 
in table[TTl and they are all in the range 0.13 to 0.23. The 
definition used here is the standard 



a = 



2E, 



NX 



ki (ki 



1 



(3a) 



(3b) 



where N is the number of nodes in the network, -Enn 
is the number of links between nearest neighbors of node 
i and fc, is the degree of node i [14|. A network embed- 
ded in two-dimensional Euclidean space with no crossing 
links is called a planar network, and is described by West 
17j . Generating a planar randomized network for com- 
paring the cluster coefficients is very hard since no links 
can cross and the rewiring therefore must be local. How- 
ever the clustering can be compared with the Delaunay 
network [l8| for the same spatial layout of nodes. For 
a given spatial node configuration and degree distribu- 
tion the Delaunay network gives the maximum possible 
clustering coefficient. The clustering coefficient for a De- 
launay network made from nodes randomly distributed 
in the plane and with a number of nodes comparable 
to our samples is 0.44. Delaunay networks are closely 
linked to Voronoi diagrams, and both are described bel- 
low. The cluster coefficient for a non-planar random net- 
work, where the links can cross, having the same number 
of nodes and links and the same degree distribution is 
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in the order of 0.001. The ridge networks have a much 
higher clustering than the non-planar networks. This is 
expected because any node in a planar network has a 
low chance of being linked to a far away node. This will 
generally increase the local clustering 17]. On the other 
hand the clustering is significantly lower than in the De- 
launay case. This indicates that the ridge networks does 
not form highly interconnected cliques. 



GEOMETRICAL PROPERTIES 

Various geometrical properties of crumpled thin sheets 
have been investigated in the past [6j, llEj, 20] . Here we 
discuss the size distribution of facets formed by the ridges 
and of the Voronoi sections formed by the location of the 
nodes. The angular distribution of the ridges and the 
3-cone structures are also investigated. 



Facets 

The nodes and links of the network form facets (also 
called domains) of various sizes and shapes. A facet is 
defined as an area of the crumpled paper confined by a 
closed loop of ridges that is simply connected, meaning 
that it contains no internal facets. The nodes bordering 
the facets are the corners or vertexes of the facet. The 
distribution of facet areas and number of vertexes for 
each sample have been calculated. The vertex distribu- 
tion for all the samples was averaged giving each sample 
equal weight. The number of facets with 3, 4, 5 and 6 
vertexes was 46%, 28%, 15% and 8% respectively, and 
the number of facets with more than 6 vertexes was 4%. 
The maximum number of vertexes was 14. In figure [7| 
the distribution of the facet vertex number can be seen, 
the data is fitted with a log-normal function 



p(a) 



1 



-(i«(Q)- M ) 2 /(2a) 2 



2naa 



(4) 



where a is the vertex number, fj, is the logarithm of 
the average number of vertexes per facet and a is the 
standard deviation. The best fit was achieved with a = 
0.42 and n = 1.13. 

The areas of the facets have also been investigated. 
The binned distribution of areas was normalized by the 
maximum area for each sample, and the average over 
all samples calculated. The resulting distribution can be 
seen in figure [8] together with a log- normal fit. The best 
fit parameters were a = 1.17 and fi — 2.16 in arbitrary 
units. 



o Vertex distribution 
Log-normal fit 




FIG. 7: (Color On-line) Plot of the average vertex distribu- 
tion, p(v), as a function of the number of vertexes v for the 
facets formed by the ridges. The data is fitted to a log-normal 
function with a — 0.42 and fi = 1.13. 
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FIG. 8: (Color On-line) Plot of the average area distribution, 
p(a), as a function of the facet area a for the facets formed by 
the ridges and the Voronoi regions. Both data are fitted by a 
log-normal function with a — 1.17 and fj, = 2.61 for the facets 
and a — 0.74 and /j, — 2.73 for the Voronoi regions, both in 
the same arbitrary units of area. 



Voronoi networks 

Given a set of nodes in space (or the plane) the Voronoi 
diagram (lij is a sectioning into areas around each node 
where each section contains all the points that are closest 
to the node in its interior. This partitions space (plane) 
into sections filling the whole space (plane). The Delau- 
nay network is the network where each node is linked to 
all the other nodes that it shares a Voronoi section border 
with. A visualization of this is given in figure [5J where 
the Voronoi diagram for four of the samples are plotted. 
The coloring of a given Voronoi section reflects the size 
of the section. Smaller sections have a brighter shade 
and larger sections have a darker shade. It can be seen 
that the sections are grouped according to size, making 
regions of the whole diagram that contains mainly large 
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FIG. 9: (Color On-line) Plot of the Voronoi diagram for four 
samples. The color coding of the various sections represents 
the area of the section. Brighter areas are smaller. There is 
a clear trend for sections of small (large) size to group with 
other small (large) sized sections. 



or small sections. The distribution of the areas of the 
various Voronoi sections have been calculated and fitted 
with a log-normal function. As in the facet case each 
sample has been normalized by its maximum area. The 
distribution follows the same general shape as the facet 
distribution, and they can both be seen in figure [SI 



Angular distribution 

The angular distribution of the ridges relative to the 
border of the sample has been studied in order to de- 
tect any preferred ridge direction or ordering among the 
ridges with regard to direction. No such preferred direc- 
tion or ordering was found, and the distribution of ridge 
angles was reasonably uniform, both for each sample and 
for the average. A plot of the binned ridge angle distri- 
bution can be seen in figure I10I 

The distribution of angles between ridges in a three 
ridge cone (a node where three ridges meet, and hence 
form a cone like structure) have earlier been investigated 
both analytically and experimentally d, 0, [2(| . It has 
been reported that there is indications of preferred open- 
ing angles for such cones in the regions about 20°, 60° 
and 110°, although all acquired distributions have been 
broad. All k = 3 nodes have been investigated and the 
distribution of the ridge separation angles shows a broad 
distribution of angles with a maximum in the range be- 
tween 100° and 150°. There are no significant peaks 
in the distribution and this indicates a random order- 
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FIG. 10: (Color On-line) a) Plot of the distribution of ridge 
angles, p((j>), as a function of angle (j> for all the samples and 
the average (elevated for clarity). No trend is visible in the 
plot, b) Distribution of separation angles, p(<f>), as a function 
of opening angel in 3 ridge nodes. The average is elevated 
for clarity. 



ing. However 32 % of all the angles lies in the interval 
between 90° and 150°. This suggests that the ridges 
tend to span out trying to separate themselves from each 
other. Recall that 120° is the angle at which they are 
evenly separated. A plot of the distribution for all the 
samples and their average can be seen in figure [TOl 



ROUGHNESS 



The roughness of crumpled paper surfaces has been 
investigated before 0, H, E3] • These investigations have 
reported self-affine behavior, this means that the surface 
is statistically characterized by 



h(x) = \- H h(\x) 1 



(5) 



where h(x) is the height of the profile at position x, A 
is a rescaling factor and H is the Hurst exponent. We 
have investigated the one-dimensional profiles produced 
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CONCLUSION 
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FIG. 11: (Color On-line) a) Results from the AWC analysis 
showing the average wavelet coefficients, Wt,{a), as a function 
of scale a, the scale represents measured points, b) Results 
from the PSD analysis showing the power spectrum density, 
P(k), as a function of spatial frequency k. In both a) and b) 
the data sets are vertically shifted for clarity. 



by the profilometer using the average wavelet coefficient 
(AWC) method [2l|, the power spectrum density (PSD) 
method [13] and the bridge method [23| . The results from 
all the methods indicates that the crumpled paper forms 
a self-affine surface. Earlier works have reported a small 
scale region with a Hurst exponent H5 ~ 1.0 and a large 
scale region with H L ~ 0.7 [EEl and H L - 0.8 [4]. Our 
results follows the same trend in that there is a cross over 
scale between two scaling regimes. However we found the 
small scale exponent to be Hs = 1.25 ± 0.05, indicating 
that the surface is asymptotically non-flat at these scales. 
Unfortunately the data did not give a robust value for 
Hi because the sample size was too small compared to 
the crossover scale. The data did however indicate that 
Hl < 1.0 and in the range reported above. In figure [TTl 
results from the PSD and AWC methods can be seen. 



The main points reported above are that the tail of 
the ridge length distribution is found to be well repro- 
duced by a power-law distribution, and that the short 
ridges follows a log-normal distribution as reported ear- 
lier. The degree distribution has been shown not to have 
a power-law tail, but rather an exponential decay, and 
the networks have been found to be disassortative. The 
facet area distribution, the corresponding Voronoi dia- 
gram area distribution and the Delaunay vertex distribu- 
tion have all been found to fit log-normal distributions. 
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